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Abstract 

A relativistic quark model of mesons formulated within the formalism of Fokker-type action 
integrals is proposed, in which an interquark interaction is mediated by scalar-vector superposition 
of higher derivative fields. In the non-relativistic limit the model describes a two-particle system 
with the linear potential. In order to analyze the model in the essentially relativistic domain the 
■ perturbed circular orbit approximation and certain principle of selection of physically meaningful 

r~| . solutions are applied which permit one to perform the canonical quantization of the model. It is 

shown that the model reproduces well specific features of the light meson spectroscopy. 
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I. INTRODUCTION 



It is known that spectra of heavy mesons (contaning c and b quarks) are described well 
by means of potential models with the non-relativistic Cornell potential u{r) = UQ — a/r + ar 
and various quasi- relativistic corrections of scalar- vector type iH^. The potential is QCD- 
motivated: its Coulomb part is a non-relativistic limit of the one-gluon exchange interaction 
while the linear part comes from the Wilson loop. The latter is also related to a string 
conception of hadrons 5|-[3]. Constants Mq, a and a vary from one model to another. In 
particular, the string tension parameter a is frequently used as an adjustable parameter from 
the range a = 0.15^0.3 GeV^ il 4] although the most conventional value a = 0.18^0.2 GeV^ 
is substantiated by QCD simulations on the lattice 8|. 

Mass spectra of light mesons (consisting of u, d and s quarks) possess characteristic 

nnn 

features which can be summarized roughly in the following idealized picture p, 191-11 Ij : 

1. Meson states are clustered in the family of straight lines in the (M^,j)-plane known 
as Regge trajectories. 

2. Regge trajectories are parallel; slope parameter a is an universal quantity, a = 1.15 
1.2 GeV^ 

3. As states of quark-antiquark system mesons can be classified non-relativistically, by i 
and Ur (the orbital and radial quantum numbers) as well as by s and j (the total spin 
and angular momentum). 

4. Spectrum is £s-degenerated, i.e., masses are distinguished by i (not by j or s) and n^. 

5. States of different i and possess an accidental degeneracy which causes a tower 
structure of the spectrum. 

Items 1-4 imply that in the (M^,£)-plane meson states form into strait lines too: the 
principal (n^ = 0) Regge trajectory built up of the set of degenerated singlet (s = 0) and 
triplet (s = 1) states, and the family of daughter trajectories {ur = 1,2,...). Hence energy 
levels of q-q system can be described by a formula: 

^ a{i + xrir + C), (1) 
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where the intercept constant ( depends on a flavor content of mesons f« 1/2 for (vr-p)- 
family of mesons; it grows together with quark masses). Finally, the accident degeneracy 
(the item 5) constraints the constant coefficient x determining the spacing of daughter 



12| 



trajectories to an integer or a rational number 

Light mesons are essentially relativistic two-quark systems, and considerable amount of 
various relativistic models has been invented for their description. The most elegant and 
listorically important among them are the simple relativistic oscillator (and its variations) 
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15| and string models |5-7|. They lead exactly or asymptotically (at large i) to the 



formula ([T]) with x = 2. Besides, the string models tie the slope and the string tension 
together, a = ka, with the slope coefficient = 27r ^ 6.3, so that, the value a = 0.18 GeV^ 
is preferable. 

Farther relativistic potential models are based on various relativistic generalizations of 



Schrodinger equations with a confinin. 



one- 
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18| and two-particle Q, ll[ 



Cornell or more complicated}__potential such as 



23| Dirac equations etc |2J, |25|. These models 



incorporate the description of heavy and light hadrons and, in most, reveal asymptotically 
linear Regge trajectories ([1]) of the slope a = ka with the slope coefficient k = 4-^8 and with 
the daughter spacing coefficient x = 2. With this kind of degeneracy, i.e., of (£+2nr)-type, 
however, a certain number of states falls out the description 0, [lO[. The value x = 1 is 
more adequate to experimental data. In particular, it follows (in the limit i ^ 1) from the 
mass formulae derived by means of the Dirac- type equation model {25] and selection rules 
superinduced by hands, and used for a description of vr-, p- 26| and K-trajectories {2^. 

In the present paper we consider the relativistic potential model of mesons which reveals 
asymptotically linear Regge trajectories with native (£+nr)-degeneracy. A classical proto- 
type of the model was formulated independently by Rivacoba 28(] and Weiss 29|| by means 



of the Fokker-type action integral 
field theory 
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3l| related, in turns, to a higher-derivative gauge 



33[. Namely, the interaction between particles is described in terms of a 



time-symmetric Green function of a fourth-order field equation. 

Hamiltonization and quantization of Fokker-type systems is rather challenging problem 
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37|. 
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le Hamiltonian description 



36| which, in most, are not 



in view of a time-nonlocal character of the interaction 
in this case can be built by means of approximated methods 
appropriate for strongly coupled systems. 

For particular time-symmetric Fokker-type systems one can invent naturally time- 



asymmetric counterparts in which a time-nonlocahty is removed [38[. The Rivacoba- Weiss 
model is the case. For this but time-asymmetric model an exact Hamiltonian formulation 
see |38l| for general formalism) and the corresponding quantum description was elaborated 
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39[. In despite of an admired degeneracy (i.e., with x = 1), the slope of asymptotic 
Regge trajectories turned out to be overestimated, with the coefficient k = 3\/6 ~ 10.4. The 
reason perhaps resides in the fact that the vector character of interaction brought into the 
model from the underlying gauge theory is not quite suited to an actual nature of a confining 
interaction in hadrons. Unfortunately, the Fokker-type model of scalar confinement without 
time-nonlocality is unknown. 

Recently, a quantization method of two-particle Fokker-type systems in an almost- 
circular-orbit (AGO) approximation has been proposed by the author j4^. The method is 
appropriate for strongly coupled systems. Here it is applied to a quantization of the time- 
symmetric Rivacoba- Weiss model. Moreover, the analogue of the Rivacoba- Weiss model 
with scalar confining interaction is built, and the scalar- vector superposition model is con- 
sidered. It is studied an asymptotic behavior of the Regge trajectories, from which the slope 
and daughter spacing coefficients are found and compared with data from experiment and 
other potential models. 

II. VARIOUS FORMULATIONS OF FOKKER-TYPE ACTION INTEGRAL WITH 
A VECTOR LINEAR CONFINEMENT. 

We start with the manifestly covariant two-particle action 

2 r r- 

I = /free + /int whcrC /free = / dx^ ^ ^ (2) 

a=l ^ 



and /int is the Fokker action integral 30|, |3l| describing an interaction. For the arbitrary 
interaction of a vector type we have: 

//nt = -j I dr2 Xi-X2 G^xl^)- (3) 

In eqs. ([2]) and ([3]) is a rest mass of ath particle (a=l,2); a;^(Ta) (/i=0, 3) are covariant 
coordinates of a world line of ath particle parameterized by an arbitrary evolution param- 
eter Ta, X^(ra) = dx'^/dTa] = ^li'Tl) " 2:2(^2); ^2 = Vf^iy^ 12^12', ^hc fuUCtloU ^(x^s) 

is usually proportional to a symmetric Green function of an appropriate field equation, 



or it may be chosen phenomenologically. We use the time-hke Minkowski metrics, i.e., 
II Vfii^ II = diag(+, — , — , — ), and put the hght speed to be unit, c = 1. 

If one chooses ^(0:^2) ^(^12) where 5(a;^) is the symmetric Green function of the 
d'Alembert equation = 4:Tt6{x), one arrives at the Wheeler- Feynman electrodynamics 



4l|. 



Let us consider the Fokker-type action proposed by Weiss [29|. It corresponds to the 
choice (j(a;^) oc 6(a;^) in ([3]) (where Q{x) is the Heaviside step function) with some coefficient 
of proportionality which we specify here as follows: 



iae(x2). 



a > 0. 



(4) 



In the non-relativistic limit the Weiss action leads 42| to the interaction potential: 

POO POO 

U{r)= MG{^f-r^) = -a dt9 = a(r - 00) (5) 



which corresponds to a linear confinement up to unessential infinite constant. 
As it is shown in 



32 



33| the Weiss action principle is related to the higher-derivative 



theory of the vector field proposed by Kiskis 43|] and to its later non-Abelian version 44. l45| 
In particular, the function B(x^) is a symmetric fundamental solution of the equation: 



D^e{x^) = 16n6{x). 



(6) 



The Fourier trans 
gluon propagator 



brm of this solution oc coincides with the infrared asymptotics of 



44|. 



An infinite constant in r.-h.s. of ([5]) indicates that the Fokker action integral ([3]) with 
the Green function (jlj) is not well posed from the mathematical viewpoint. A formal causal 
structure of the interaction is that as if each point (say, Xa) of a world line of one particle 
is related to infinite segments of another word line lying inside the light cone with the 
center Xa, and the contribution of these segments in the action is infinite; see Fig. 1(a). 
Physically it is not crucial since a variation of the action turns 6(a;^) into its derivative 

G'(x^) = S{x^), and Euler-Lagrange equations relate points of particle world lines along 
generatrices of light cones only; see Fig. 1(b). Nevertheless, integrals of motion such as the 
energy and the angular momentum turns out divergent. In order to avoid this difficulty one 



can reformulate the Fokker action ([3]), via the integration by parts 46| : 



'int 



/oo /»oo 
/ dri dr2 Xi-X2 6(x 
-oo J —oo 



12J 



00 ^ —00 

00 POO 



dridr2 (xi2-xi)(xi2-X2) (^(a; 



12) 



(7) 



00 ^ —00 



^©(^12)^12 



ri=oo 
ri = — 00 



r2=oo 
r2 = — 00 



The last divergent term does not contribute in the equations of motion, and we arrive at the 
equivalent formulation of the problem proposed earlier by Rivacoba [28]. The Fokker- type 
integral ([7]) itself describes an interaction with the causal structure of Fig. 1(b), as in the 
Wheeler-Feynman electrodynamics, and leads to finite integrals of motions. 




Figure 1. Interaction causal structure of various Fokker- type action integrals (specified in the text). 
Solid curves depict world lines of particles 1 and 2. Arrows and thin lines depict generatrices and 
inwards (a) or outwards (c) of light cones where points of particle world lines are related. 

One can propose third equivalent formulation of the problem which is most convenient 
for our purpose. Using the equality 9(a;^) = 1 — 0(— a;^) one obtains: 

;>oo ;>oo 

4\^ = ^/ / dridr2Xi-X2 [l-e(-xy] 

J —00 J —00 

= -- / / dridr2 a;i-a;2 6(-a;i2) 
J —00 J —00 



(8) 



-CXD ^ —00 

Tl=00 



Tl= — 00 



r2=oo 
r2 = — 00 



The interaction causal structure of the integral ([H]) is shown in Fig. 1(c). The integrals of 
motion are the same as in the Rivacoba version of the model. This version (i.e., (jS])) of the 
Weiss action ([3]), (jl]) is equahy substantiated by the Kiskis field theory since the function 
— 6(— x^) satisfies the equation ([6]) as well. 
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III. AGO APPROXIMATION IN THE FOKKER-TYPE DYNAMICS. 



The construction of the Hamiltonian description of Fokker-type systems, as a step towards 
quantization, is a rather difficult task which can be reahzed within a certain perturbation 



scheme. A commonly used quasi-relativistic approximation scheme (see, for example, 34| ) 
works well if relativistic effects are weak. But light mesons, as two-quark systems, are 
essentially relativistic, and they need another approach. 

Here it is used the almost-circular-orbit (AGO) approximation scheme developed in the 
previous work of the author [401]. The scheme is based on the fact that all known in literature 
two-particle Fokker-type systems with attractive (in some meaning) interaction posses exact 
solution of the shape of concentric planar circular particle orbits of radii Ra{^) dependent 



on an angular velocity Q; see 
Fokker-type system: 



^7H49|. In |40| this is proven for a general two-particle 



2 r 
a=l ^ 

+ j j dhdt2^{tiM,Xl{ti),X2{t2),Xi{ti),X2{t2)) (9) 



which is invariant under the Aristotle group (including time and space translations and 
inversions, and space rotations), at least. Manifestly covariant Fokker-type systems 
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4l| which by construction are Poincare-invariant (and the more Aristotle-invariant) 
as well as possessing reparametrization invariance can be reduced to the form ^ by means 
of the choice of the evolution parameter = ta = 2;°; then the particle positions are 
Xaita) = {^aita)} (i = 1,2,3). The manifestly covariant Rivacoba- Weiss system ([2]), ([7]) 
or ([8])) does possesses exact circular orbit solutions even in a strongly relativistic domain 



Thus a set of these solutions can serve as a zero-order approximation in a perturbative 



treatment of the Fokker-type dynamics. 

The invariance of the action ([9]) with respect to time translations and space rotations 



leads to an existence of the energy an the angular momentum integrals of motion |50| : 

^ = t|--al:-i}(i" + A")+#dMt.{A-A|t. (10) 

9 

dXa 



a=l 
2 



d 

a=l 



ijjf dt, dt2 { + 0.2) X ^ + X A _ X A I (11) 



where 



/oo rCG nn rt\ poo poo nt2 

dt2$, A2= / dti*, K = - . 

oo J —oo JJ J —oo J t2 J ti J —oo 

On the circular orbits these integrals are functions of the angular velocity: -E'(o)(^^) and 
J(o)(r2), J(o)||r2 so that we can get J(o)(^^) where J(o) = | J(o)| and Vt = \ fl\. 

Let us transit to a non-inertial reference frame which is uniformly rotating with the 
angular velocity O. This can be done via the change of variables Xa(ta) — )• Za(ta): Xaita) = 
S(ta)za(ta) wherc S{t) = exptQ G S0(3) and the skew-symmetric matrics Q is dual to the 
vector n. Within this reference frame a circular motion of particles is described by static 
vectors Ra such that R2 tJ- Ri- Then small perturbations of circular orbits are characterized 
by deviation vectors Paiia) = Za{ta) — Ra- 

Expanding the action ([9]) in powers of pa yields in the lowest non-trivial order the 
quadratic form: 

^^"^ = l5^//dtdtV^'(i)^fcKi-0p'(0, (12) 

kl 

where the kernel matrics D(t — t') = \\Dki(t — is invariant under time translations and 
reversion: D^(t' — t) = D{t — t') (here the multi-indeces k^l = {a,i), {b,j) has been used). 
Corresponding equations of motion form a time-nonlocal linear homogeneous system: 

J2fdt'Dki{t-t')p'{t') = 0, (13) 
I 

which possesses a certain fundamental set of solutions. Among them the exponential solu- 
tions p^(t) = e^(a;)e~"^* are of interest. Substituting them into the system f|T3|) yields the 
set of algebraic equations: 

J2Dkiioj)e\u)=0, (14) 

which amounts the eigenvalue-eigenvector problem for the polarization vector e^(a;) and 
the frequency u. The latter is determined by means the secular equation det D(a;) = 
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in terms of the dynamical matrix D{uj) = Jrft D(t)e''^*. In view of time-nonlocahty of 
the problem f|T3l) the matrix entries Dki{uj) are, in general, non-polynomial functions of 
u, and the set of solutions of the secular equations may be infinite. Due to symmetric 
properties of the dynamical matrix this set consists of duplets if € K or quadruplets 
{±UJa, ic^a, a = 1, 2, . . . } if Im 7^ 0. In the latter case the corresponding solution is 
unbounded and cannot be described correctly within AGO approximation (where p'^ must be 
small). Thus among all eigenfrequencies we select real ones only and arrive at the following 
solutions of the system ( IT3|) : 

p\t) = J2 {Aae^M e-'--' + AlZ{u^) e'-Stj , (15) 

a 

where complex amplitudes Aa of oscillations (modes) parameterize the phase space of the 
system. Only one mode corresponding to mutual radial particle oscillations with the 
frequency Ur is physically meaningful. Other modes are either kinematic ones which can 
be reduced via redefinition of zero-order circular orbits, or non-physical ones which reveal 
physically unacceptable behavior of particles and arose as a mathematical artefact of the 
theory (as in the Lorentz-Dirac equation, for example). All such modes should be discarded. 
After this is done and the polarization vectors e^{(jJa) in ( IT5|) are appropriately normalized, 
the angular momentum and the energy of the system take the form: 

j = Mm. (16) 

E = Ei^o){9) + E^2){^,Ar) (17) 
where Et^i^iVt, Ar) = ujr{yt)\Ar\^ . (18) 

Other integrals of motion following from the Poncare-invariance of the system vanish; they 
are the total momentum, P = and the center-of-mass integral (boost), K = Q. Thus the 
AGO approximation brings the system into the center-of-mass reference frame. 

In order to construct the center-of-mass canonical description of the system one should, 
first of all, to invert the relation (|T6l) with respect to Vt = Q{J). This permits us to obtain 
the center-of-mass Hamiltonian which is nothing but the total mass of the system: 

M = M(o)(J) + M(2)(J, 1^1) = {E^om + cOrmArn^^^^j^ . (19) 

It is understood as a function of J = \J\ where components Jj {i = 1,2,3) of the intrinsic 
angular momentum J of the system satisfy the Poisson bracket relations (PBR): 

{ Jj, Jj} = eij'^Jk, (20) 
9 



and of the amplitude of interparticle radial oscillations Ar satisfying the PBR: 



{Ar, a;} = -i , {Ar, Ar} = {A^, A*,} = 0. (21) 

In order to transit to an arbitrary reference frame one must introduce canonical variables 
characterizing the state of the system as a whole, for example, the total momentum P and 
the canonically conjugated CM position variable Q. Then a complete Hamiltonian descrip- 
tion of the system, i.e., ten canonical generators of the Poincare group, are determined in 
;erms of M, J, P and Q via the Bakamjian-Thomas (BT) model or equivalent constructions 
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52| . The quantization of BT model is well elaborated 53l. |54 1 . 
In present work we are interested mainly in the spectrum of the mass operator M. It can 
be obtained directly from (fT9|) by means of the following substitution: 

Ar ^ Ar, A >■ AL] 



J ^ VJ2^ £ = 0,1,...; (22) 

\Ar\^ ^ \{Aral + alAr) ^ TLr + \, n^=0,l,... (23) 

Here the condition ra^ ^ £ is implied, due to a perturbation procedure. 

IV. RI VAC OB A- WEISS MODEL IN AGO APPROXIMATION. 

Let us consider a circular-orbit solution of the Rivacoba- Weiss model. Using the action 
(E]), dH]) for a system of two equal particles of the mass = m (a = 1, 2) and following 
the general methodology proposed in 49| or [4^, one states a relation between the angular 



velocity of a motion of particles along circular orbits and the radius R of these orbits. 
It is convenient, instead of R, to handle with particle velocities Va = v = Rfl. Then the 
relation between Q and v can be determined implicitly, or parametrically, via an auxiliary 
angle 0. It is related with the velocity v by means of the equality: 

/((/)) = - 4v^ cos\(j)/2) = (0 < i; < 1) =^ 

9 02 6/2 , , 

2(1 + cos 0) cos(0/2)' ^ ' 



In turns, we have for f2 

(25) 



m 
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where /'(</>) = d/((/)) /d</), F = (1 — t>^)~^/^ and the superscript "(v)" refers to the vector 
interaction. Let us note that v E (0,1), R E (0, oo) and Q G (oo, 0) if G (0, 0i) where the 
angle 0i/2 = xi = 0.235 vr is a positive solution of the transcendental equation x = cosx- 

The integrals of (circular) motion M(o) and J = J(o) are convenient to write down as 
follows: 



a 



1 + - 



f(0) 



= (26) 



- m^/') - (27) 



a 



They grow as M(o) G (2m, oo) and J G (0, oo) if G (0, 2xi). 

In order to study the system in AGO approximation we need to construct the reduced 
2x2 dynamical matrix [40] and then to calculate the frequency lOr of radial oscillations 
or, equivalently, the fraction A = uJt/VL, as a function of either J, v or (which is most 
convenient) 0. This is done in the Appendix A. 

Here we are interested of an asymptotic expression for the total mass f|T9|) squared at 
J — 7- oo. Within the perturbation procedure the inequality M(2) <^ M(o) is implied. Taking 
this into account one obtains: 

^ Mj) + 2M(o)M(2) = Mj) + 2M(o)u;r.|A,|2 



If the following limits 



, U.^^.™|f. (29) 

J-i-OO J (p-Kpl Jj[(p) 

X = 2 hm -^A = 2 hm |4Sa(0) (30) 



exist and are finite, the asymptotic value for the total mass squared (!28l) takes the form 

~ /ca{J + at J^oo (31) 

and, upon quantization (l20l) . recovers the Regge trajectories (P) with cr = ka. In the present 
case of the vector confinement model 

fcW = 8xi(l + sinxi) ^ 9.896, x^^^ = 1. (32) 
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An asymptotic value of the daughter spacing coefficient x'^^) = 1 matches well for a de- 
scription of the tower structure of meson spectra (see item 5 in Sec. 1). But the slope 
coefficient k^^^ ^ 9.896 exceeds conventional values /c = 4 8 which occur in various po- 
tential models. A plausible reason of this disagreement in that the purely vector nature 
of interaction in the model does not correspond to the actual relativistic structure of the 



confinement interaction which is commonly opined as of scalar 
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m 



22| or scalar-vector 



19| type. 



In order to confirm or challenge this assumption we construct in the next section the 
scalar analogue of the Rivacoba- Weiss confinement model. 

V. THE FOKKER-TYPE ACTION INTEGRAL WITH SCALAR LINEAR CON- 
FINEMENT. 

Let the two-particle action to include the free-particle terms ([2]) and the Fokker-action 



integral of a scalar type 



33|: 



r(s) 



J j dri dr2 \f^i\f^G{x\2). 



(33) 



If the function G(x^) is chosen in the form (jl]) it is expected that the action ([2]), fl55]) 
describes the scalar confinement interaction. Indeed, the action fl33p can be derived from 
the higher-derivative theory of scalar field jsS ]. 

In this case however one encounters even more significant divergences as in the vector- 
type model since not only the action itself and integrals of motion but also the equations of 
motion are ill-posed. Fortunately, the remedy to set the scalar model properly is the same: 
one replaces the function (jl]) by 



G{x^) = ia0(-a;2), a > 0, 



(34) 



which is analogous to the transition from the action (jl]) to (|H]) in the case of Weiss 
model. The replacement of the function (jlj) by (I34p in the action (13 3 p may also be treated 
as a renormalization of particle rest masses: 



a = 1, 2, a = 3 — 1, 



(35) 



where mo a is an infinite bar mass of ath particle and nia is finite. 
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A subsequent consideration of the scalar model is similar to one in the vector case. The 
system of equal rest masses is considered. Dynamical characteristics of circular orbit solution 
are parameterized by the angle (j). In particular, for the angular velocity VL one can obtain: 



a 



1 



(36) 



where f'{(f)), v and F as functions of are defined in Sec. 4. In contrast to the vector case, 
here f^\4>) — )■ if — 0o 7^ 0i where 0o/2 = Xo is a positive solution of the transcendental 
equation: 

3x' cos X + sin x - cos^ X = 0, X e [0, 7r/4] . (37) 
The latter by means of the substitution x = 7r/2 — S'?/' can be reduced to the form: 

V^ = ^-|^, ^G[vr/12,7r/6] (38) 
6 ocos-^ 

which is convenient to iterate the numerical solution: Xo = O.lSlvr < Xi- It is surprisingly 
that particle velocity v — 0.535 < 1 at x — ^ Xo while orbit radus i? — t- oo. This distinguishes 
the scalar model from the vector one in which f — )■ 1 at i? — t- oo. For the integrals of (circular) 
motion we have: 

nM(o) 20^(1 -.^)^ _ ,) 

^ = (1-^^)0-/^(0)- (40) 
a 

It is easy to verify that M(o) G (2m, oo) and J G (0, oo) if G (0, 2xo)- 

Using the functions fl5^ . f HU]) in eqs. fl2^ . f l5U]) and taking limits at — )■ 0o (instead of 
— )■ 0i) yields the slope and daughter spacing coefficients: 

A;(^) = 2.716, x^*^) = 1.902. (41) 

The latter is close to 2, as in the oscillator-like and some string relativistic models of mesons 
[7. [131 15|. The accidental degeneracy and thus the tower structure of the mass spectrum 
is recovered approximately. Again, the slope coefficient is not appropriate (similarly to the 
vector model), but it is considerably less than the conventional values 4-^8. 

The difference between the vector and scalar models suggests that general features of the 
light meson spectroscopy may be recovered (at least, asymptotically) within the Fokker-type 
model with a scalar-vector confining interaction. 
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VI. THE FOKKER-TYPE ACTION INTEGRAL WITH A SCALAR- VECTOR 
CONFINING SUPERPOSITION. 



The Fokker-type system of two particles bound via superposition of scalar and vector 
confining interactions is naturally defined by means of the action ([5]) with 



'mt 



r(v) 



(42) 



where I^^^ and I^^^ are defined in eqs. (ES]), f l51|) and ([H]), respectively, while ^ G [0, 1] is a 
mixing parameter. 

All the functions f^\(j)), and fj^\4>) determining the dynamics and integrals of 

circular motion of this model are superpositions of the functions fj25|) - fl271) and fl36l) . fl39|) . 

/(«)(0) = (1 -0/(^^0) (43) 
Then the function [M^q}^{J)]'^ which is a classical analogue of the principal Regge trajectory, 
can be presented in the parametric form: 



r(v) 



M, 



(0) 



aJ 



n 2 



2 ' 



[0,1]; 



(44) 
(45) 



it is shown in Fig. 2. The maximal angle 0^/2 = Xi is the smallest positive root of the 
equation f^\2x) = 0. It grows monotonically over the segment ^ [XcXi] if ^ ^ [O5 1/2], 
and = Xi if C ^ [1/2? 0]. Similarly, the maximal speed of particles (at -R — 00 when 
M(o) — !■ 00 and J — t- 00) grows monotonically, v G [0.535,1] if ^ G [0,1/2], and t> = 1 if 
^G [1/2,1]. 

The slope and daughter spacing coefficients can be calculated similarly to the previous 
cases, i.e., using eqs. (|29|) and f l30|) with the limiting angle (p^ (instead of (pi). One can proof 
that the following equality holds: 



/f (0) 



lim 



[0,1]. 



Thus the formula (I5U1) for the daughter spacing coefficient simplifies: 

= lim X^^\(j)), e G [0,1]. 



(46) 



(47) 
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lb 20 30 40 50 

Figure 2. Classical Regge trajectories for different values of the mixing parameter ^. 

The function A*^^^(0) is determined numerically from the secular equation detP^^^(A) = 
for the matrix flA.171) : see Appendix where the graph of X^^^cp) is presented in Fig. 5. 

Both the slope and daughter spacing coefficients are functions of the mixing parameter. 
In particular, 

= ^k^^\ x(«) = 1, e e [1/2, 1]. (48) 

A behavior of these functions on the whole segment C, € [0, 1] is presented in Fig. 3. Grid 
lines on the graphs take values of ^, k and x into a mutual accordance for particular cases 
^ = 1/2,1 and k = 4,27r,8. 




Figure 3. Slope (left graph) and daughter spacing (right graph) coefficients vs mixng parameter in 
the scalar- vector model. Grid lines connect values ^ with k and x at ^ = 1/2, 1 and k = 4, 27r, 8. 

It is seen from these graphs that the slope coefficient k^^^ is a monotonically increasing 
function of the mixing parameter ^: k^^^ E [2.716,9.896] if ^ [0? !]• This segment includes 
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conventional values of = 4 8 which occur in non-relativistic and relativistic potential 
models. 

Degeneracy properties of the system with scalar-dominating confinement interaction (i.e., 
at ^ < 1/2) differ crucially from those of ^ > 1/2 case. In particular, the vector-dominating 
model possesses the asymptotic accidental degeneracy of {i+nr)-tjpe. Since one can provide 
in this case an arbitrary value for k from the segment k^^^ G [4.948,9.896], the vector- 
dominating model may be compared to variety of non-relativistic potential models and 
string model. 

For the scalar-dominating model the lower conventional bound k = 4 for the slope is 
achieved at the mixing ^ ^ 0.37 which, in turns, leads to the daughter spacing x ^ 3/2. 
The accidental degeneracy is present but somewhat hidden in this case. 

Upon quantization of the model the mass squared spectrum is calculated by means of the 
quantization rules ( 122|) -( 123|) used in the classical expression f l28|l . Practically, one substitutes 
J = i + ^ in l.-h.s. of ( l45l) and solves this equation for angles (pi {£ = 0, 1...) which, in turns, 
are used as arguments of the functions (14^ . fn{4>) and A(0) in r.-h.s. of (l28l) . 

"1 £ 

rir — 1 2 3 4 




2 4 6 8 10 12 14 16 18 



Figure 4. Quantum Regge trajectories for a = 0.18 GeV^, m = 0.15 GeV and = 0.63. The 
asymptotic slope a = 27ra = 1.15 GeV^, the daughter spacing x = 1. 

Let us note that classical Regge trajectories (Hljl . (and Fig. 2) start from J = 
corresponding to = 0. In the quantum case the bottom value for the dimensionless quantity 
j = Ja/rn^ (in l.-h.s. of ( l45ll ) corresponding to s-states (i.e., £ = 0) is jo = ^a/fn'^ > 0, hence 
00 > 0. For example, taking a = 0.18 GeV^ and m = 0.15 GeV (the constituent mass of 
light quarks) yields jo ~ 4. In this case a notably curved bottom segment of classical Regge 
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trajectories which is present in classical case (see Fig. 2) disappears from the quantum 
principal trajectory which thus is closed to a straight line ([1]). Instead, daughter trajectories 
acquire an erroneous curvature in their bottom, due to an inapplicability of the quantization 
method at rir > i. This is illustrated in Fig. 4. It is seen an approximated tower structure 
of spectrum, due to the asymptotic degeneracy of (£+nr)-type. 



VII. DISCUSSION. 



In the present 



Weiss model 



28 



paper the ACO-quantization method 40[ has been applied to the Rivacoba 



29| . This model represents a Fokker-type system of two particles which 



interaction can be interpreted in terms of the classical higher- derivative theory of a vector 



gauge field 



33 



43|-|45| . The Green function oc of this field behaves as an infrared asymp- 



totics of gluon propagator 



4J] and leads in a nonrelativistic limit to the linear interaction 



potential U = ar. In the ultrarelativistic limit the model reproduces asymptotically linear 
Regge trajectories whith the slope a ~ 9.9a related rigidly to the string tension parameter 
a. The energy spectrum reveals the accidental degeneracy of (£+nr)-type which provides 
a tower structure of spectrum. Thus the quantized Rivacoba- Weiss model may serve as a 
good base for a description of light meson spectra. 

In a variety of non-, quasi- and relativistic potential models of heavy and light mesons the 
linear potential U = ar appears as a scalar (or scalar-vector) long-range part of inter-quark 
interaction. If one believes that the string tension a is a universal (i.e., flavor- free) parameter 
with conventional values in the range a = 0.15 -r- 0.3 GeV^ then the Rivacoba- Weiss model 
overestimates the slope parameter a. Since this model is purely vector, its counterpart based 
on the higher-derivative scalar field theory 55|] has been constructed. The scalar model, 
however, underestimates the slope of Regge trajectories. Finally, the family of scalar-vector 
superposition models is studied. It turned out that the slope parameter a = 1.15^ 1.2 GeV^ 
and the string tension parameter a = 0.15 ^ 0.3 GeV^ can be mutually accorded if the rate 
of the vector interaction ranges C, = 0.37 ^ 0.8. Besides, a value of the mixing parameter 
^ determines the daughter spacing parameter x. In particular, x=3/2at^ = 0.37 and 
X = 1 if ^ > 1/2, so the tower structure is also provided. 

It is worth to note that within non- and quasi-relativistic potential models the linear 
interaction is meant mostly as a scalar one. But in many relativistic models, especially 
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16 



le Dirac equation, the scalar-vector structure of a long-range interaction is 



19| . In particular, the mixture C, = 1/2, as in [lOJ], or closed values C, = 0.48-f- 



those based on 
preferable 

0.65, as in 1^, enables to reduce a spin-orbital splitting in accordance to observable values. 
The present relativistic model assures the scalar-vector structure of confining interaction 
from another viewpoint. 

In order to be appropriate for the description of both light and heavy mesons the model 
should be modified. First of all, the vector short-range interaction due to one-gluon exchange 
must be introduced. It can be done naturally via complementing the action ([2]), fj42l) by 



the Wheeler- Feynman term, i.e., by (j3]) with (^(x^ 



-a 



where a is a strong coupling 



constant. Then the model reproduces, in the non-relativistic limit, the Cornell potential. 
This modification is expected to affect some characteristics of the model in a relativistic 
regime. In particular, this may change bottom segments of Regge trajectories and decrease 
their intercept ( (see ([p) by some portion oc aa, similarly to what happens in the time- 
asymmetric model 32|, |39|]. In turns, a small intercept is appropriate for a description of 



lightest mesons A study of the model complemented with the Wheeler- Feynman term 
is beyond the scope of this work. 

Another extension of the model for a sterling meson spectroscopy is the insertion of par- 



ticle spins. One can exploit, as a guidehne, a description of spinning partic 



anti-commuting variables used in the Wheeler- Feynman electrodynamics 
tion method should be modified appropriately. 



es in terms of 



56|. A quantiza- 
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APPENDIX. CALCULATION OF AND A = LOr/^ 



It is convenient to define a dimensionless 2x2 reduced dynamical matrix 



m 



nC + ]C-S = fni(p)C + K-S 



(A.l) 
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where: 



C 



r3 + Ar -i Ar3t;2 
iArv r + Ar3 



v=<f> 



I 1 d ^ 

/'(^)d^/'(^) ' 



1 d 1 _ 



ip=(f) 



(A.2) 



(A.3) 



(A.4) 



The matrix C comes from the free-particle term of the action ([2]). The function fn{4') and 
components of other matrices K, and S depend on the interaction modeL 

For the vector (Rivacoba- Weiss) model the function fQ\4>) in defined in fl25l) . and ma- 
trices in r.h.s. of (lA.SP and (]A.4p have the form: 



/C^^) = 0, 



(v) 



"0 



^1 



l+t;2c(3+2c) 



-4t;2(l + v'^c 
(l + A^ 



2t;2s2 

(l+c)2 



-2i Xv^{l + c) 



1 
-1 



s2 



cC sS 




sS 


-cC 




- 2iA 






-sS cC 




cC 


sS 



- 2i Xv 



c(l+t;2(2+3c))C s(l+t;2(l+3c))S 
-s(l+^;2(l+3c))S (lW(c2-2s2))C 

2s(l+c)S (s2-c(l+c))C 
c(l+c)-s2)C 2scS 

{1+cfC s(l+c)S 
-s(l+c)S -s^C 



(A.5) 
(A.6) 

(A.7) 

(A.8) 



(A.9) 
(A.IO) 



where s = siny?, c = cosip, S = sin(Av9), C = cos{Xip). 

For the scalar confining interaction the function f^\(f)) is defined in (!36|l . and matrices 
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in r.h.s. of f lA.SP and (1A.4P have the form: 

r2 



(s) 



) = 2 



+ -iA(r^ + l) 

1 

l-v^{3+2c) 



1 




1 








-1 




1 







(l+c)2 
s2 



2i At;^(l + c) 



1 
-1 



^1 



-2 



- 2i A?;2 



1 -iA 
iA A2 

c(l-3t;2)-2t;2)C i(l-2t;2)sS 
-i(l-2t;2)sS (l-t;2)cC 

(l+c)C 
(l+c)C -2sS 

(l+c)2C s(l+c)S 
-s(l+c)S -s^C 




0.5 



0.8 



0.2- 



O.lir 0.27r 

, 1 , . . J — ^ 



O.Stt 



0.4ir 



0.477r (/) 



1 V 



(A.ll) 
(A. 12) 
(A.13) 
(A. 14) 



(A.15) 
(A.16) 



Figure 5. The relative frequency A = cor/^ as a function of the angle (j) and the velocity v of 
particle circular motion for different values of the mixing parameter ^. 

For the scalar- vector superposition the dimensionless dynamical matrix is constructed as 
follows: 

= (1 - ^)p(") + (A.17) 
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where ^ is the mixing parameter. The relative frequency A is then calculated as a real 
positive root of the reduced secular equation det'D(A) = 0. In general, this can be done 
numerically. 

In Fig. 5 the relative frequency A = Uy./^ as a function of the velocity v of particle 
circular motion is shown for various values of the mixing parameter ^. Let us note that 

lim = a/3 



as it must be for the nonrelativistic problem with the linear potential U = ar 40 |. 
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